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ABSTRACT 

We present results from our two year study of ground-layer turbulence as seen through the 6.5-meter Magellan 
Telescopes at Las Campanas Observatory. The experiment consists of multiple, moderate resolution, Shack- 
Hartmann wavefront sensors deployed over a large 16 arcminute field. Over the two years of the experiment, 
the ground-layer turbulence has been sampled on eleven nights in a variety of seeing and wind conditions. On 
most nights the ground-layer turbulence contributes 10% to the total visible-band seeing, although a few nights 
exhibit ground-layer contributions up to 30%. We present the ground-layer turbulence on the sampled nights as 
well as a demonstration of its strength as a function of field size. This information is combined with data from a 
MASS-DIMM seeing monitor adjacent to the Magellan Telescopes to infer the annual ground-layer contribution 
to seeing at Las Campanas. 
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1. INTRODUCTION 

The concept of ground-layer adaptive optics correction (GLAO) is to make a modest improvement in image 
quality over a large field of view (many arcmin) by compensating for the turbulence in the lowest layers of 
the atmosphere. GLAO is complementary to full-atmosphere adaptive optics correction (AO), which delivers 
diffraction- limited imaging, but only over small fields of view (< 1 arcmin). 

Recent studies of the atmosphere at astronomical sites have been conducted using a variety seeing monitors, 
and also with balloons carrying micro-thermal sensors. 1-6 A number these studies provide measuremants of 
both the total seeing and the distribution of turbulence as a function of altitiude. The new results suggest that 
typically more than half of the turbulent power in the atmosphere occurs at low altitude (< 1.0 km). This 
ground-layer is thought to arise from the interaction of moving air with local topography, which may differ from 
site to site depending on wind direction and ground contour upstream of the telescope. 

The presence of a strong, low altitude turbulence layer provides an opportunity to improve the seeing using 
an AO system with modest operating parameters. 7 Near the ground, the wind speed is low so the crossing time 
for a turbulent cell is typically on the order of ten milliseconds or more. Because the effective height of the 
ground layer is largest for the largest turbulent scales, a GLAO system favors correction of the lowest order 
modes. For atmospheric layers near the telescope, the isoplanatic angle will be large. Thus, a GLAO system 
might operate at low frequency (100 Hz) and low order (<50 modes on a 6.5-m telescope) while significantly 
improving the image quality over a large field of view (many arcmin) . Initial GLAO models predicted a gain of 
a factor of two in image size over fields of 10 arcmin or more, resulting in a factor of four gain in sensitivity for 
background limited observations. 7 ' 8 The potential of GLAO correction is evident, however most of the work to 
date has consisted of modeling based on a few atmospheric measurements. Only limited amounts of data have 
been taken at the telescope. 

Although several on-going studies reveal at least some information about ground-layer turbulence, 9, 10 none 
of the current experiments provides adequate empirical measurements of the parameters and performance of 
a GLAO system. Early laser-tomography experiments, for example, are restricted to smaller fields of view. 
Since in any event different observatory sites might be expected to exhibit different ground-layer conditions, we 
have chosen to conduct an experiment to directly measure the effect of ground-layer turbulence at the Magellan 
Telescopes. This experiment was conducted as part the Giant Magellan Telescope (GMT) project, in order to 
study the implications of GLAO for the design of GMT and its instrumentation. 



2. GLAO INSTRUMENT/EXPERIMENT DESCRIPTION 

We have developed a system of multiple Shack-Hartmann wavefront sensors to investigate the ground-layer of the 
atmosphere as seen through one of the 6.5 meter Magellan telescopes. Up to four separate wavefront sensors are 
attached to a machined plate in the telescope focal plane at the locations of stars in a few bright asterisms (both 
real clusters and random groupings) selected at convenient locations around the sky. The cameras operate at a 
frame rate of 100 Hz. The separations between stars range from a minimum of 4.5 arcminutes, set by mechanical 
interference between adjacent units, to 16 arcminutes, set by the field of view of the standard Magellan guider 
package and the size of the machined plate. The system uses off-the-shelf components to minimize cost and 
development time at the expense of sensitivity. 

The optics consist of a stock achromatic field lens and a lcnslet array which produces multiple images of 
the star directly on the detector. The A.O.A. Inc. lenslet array has a pitch of 203/im and a focal length of 5.8 
mm. The optics sample the pupil in an 11x11 square grid, resulting 90 usable subapertures which are inside the 
perimeter of the primary mirror and unobscured by the shadow of the secondary. The subapertures are 59 cm 
square at the primary mirror, which is somewhat larger than the coherence scale so there is some seeing within 
a subaperture. The original choice of subaperture size was was motivated by the idea that it might be possible 
to implement a functional low-order, wide-field GLAO system based on natural guide stars. 

The simplest algorithm for isolating the wavefront error contributed by the ground layer is just to average 
the outputs of multiple wavefront sensors distributed around the field of view. The analysis is simplified if the 
individual subapertures are accurately registered to identical pupil positions in all of the wavefront sensors. This 
requires both rotational alignment of the lenslet arrays and cameras, as well as x-y alignment of the lenslet arrays 
on the pupil. In order to accomplish the x-y alignment, the lenslet array is deliberately offset from the exact 
focus of the field lens. Each wavefront sensor is mounted on an x-y stage. Adjustment of the x-y position of the 
stage changes the positions of star images on the camera, but because the lenslet array is not quite at the focus 
of the field lens, the adjustment also changes the registration of the lenslet array on the pupil. In our wavefront 
sensors, a 1.0 arcsecond motion of the star corresponds to a 75/zm motion on the detector. The spacing between 
subapertures is 2.6 arcseconds. However the motion required to change the alignment of the subapertures on 
the pupil is much larger: a displacement of 5 arcseconds will move the subaperture position about one-half of 
the distance between subapertures. Strictly speaking, the subaperture alignment should be referenced to the 
conjugate height of ground layer. This was accomplished by mounting a special pupil mask near the center of 
the Magellan secondary mirror. Since the telescope is Gregorian the conjugate image of the secondary is located 
about 60 m above the primary. The wavefront sensor mounts include an additional motion along the z-axis 
which is used to adjust the focus of the wavefront sensor (or equivalently the pitch of the Shack-Hartmann array 
on the detector) in order to compensate for the curvature of the telescope focal plane. 

The cameras are commercial units purchased from Basler Vision Technologies. The Basler A602f cameras 
use a Micron CMOS detector (MT9V403C125STM) which has a fast readout (up to 200 full frames per second), 
but is relatively noisy (80 electron rms/pixel). The cameras have a trigger signal for camera synchronization 
and a firewire (IEEE 1394a) bus for data I/O, with 8 bits of output selectable within a 12-bit dynamic range. 
The CMOS detector has 9.9 micron pixels in a 656 by 491 format. In operation we read out a 312 by 312 
section of the chip. With this region of interest, the maximum frame rate is 210 Hz, but in practice we operate 
between 100 and 128 Hz. After readout, the data are binned 2x2 which results in sub-apertures separated by 
10 pixels, allowing for accurate centroiding and spot size measurements, and eliminating sub-aperture boundary 
crossings. Because of the high read noise, the wavefront sensors have an effective magnitude limit of 7.25 on the 
6.5 meter Clay (Magellan 2) telescope. At maximum gain, variations in bias level within a frame are significant, 
but subtraction of a bias frame leaves the data flat to a few percent. The peak efficiency of the detector is 32% 
at 550 nm, and the system has an effective wavelength of 620 nm. 

The Basler A602f camera follows an industry standard for firewire cameras (version 1.30 of the "1394 - based 
Digital Camera Specification" (DC AM) issued by the 1394 Trade Association). This allowed us to write our 
own control software on a linux-based PC using the open-source implementation of the DCAM specifications 
(libdcl394). 

Each of the 4 cameras is attached to a separate computer, with a fifth computer to provide overall synchro- 
nization. The raw data rate is substantial at 0.5 Gigabyte/min/camera, but a few minutes of data are adequate 
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Table 1. Engineering runs for the GLAO experiment at Magellan. 



to characterize the ground-layer behavior at any given time, and inexpensive hard disk drives provide adequate 
storage capacity. Exposure start times are controlled through the use of a common trigger voltage sent to all 
of the cameras. Time synchronization of the frames recorded by separate computers can be a problem because 
of local clock drift between computers and occasional dropped frames. For this reason a pattern of randomly 
spaced 1 and 2 millisecond pauses is imposed on the basic 100 Hz trigger signal, which establishes a unique time 
sync pattern which can be recovered by looking at the difference in time stamps between frames. The hardware 
setup allowed for continuous recording of data and all the frames were saved to disk for later processing. 

3. OBSERVING RUNS 

We conducted five observing runs over two years (2004 and 2005), mostly using engineering time and sharing the 
telescope most nights with other tests or science observations. The first run was used to field test the hardware 
and work on the software. The other four runs produced 11 nights of wide-field, 100 Hz wavefront sensing data. 
Typically the setup would take half an hour per camera and then once all the cameras were aligned, the wavefront 
sensing data would be recorded for approximately ten minutes. Setup involved initializing the hardware and 
then the aligning the pupil on the lenslet array in each wavefront sensor. 

Table lists the parameters of each data set. The number of cameras and sampled baselines between the 
cameras are listed in columns 4 and 5. The seeing from the Clay guider, DIMM and MASS instruments are listed 
in columns 6-8. These seeing measurements are the averages in the ten minutes prior to the observation end 
time (listed as a Julian Date in column 3). All of the seeing measurements have been corrected for airmass and 
to a reference wavelength of 0.5 [im. The wind speed and temperature are reported as the average conditions for 
ten minutes surrounding the observation. 

4. REDUCTIONS/RESULTS 

The imaging data was corrected for bias offsets, binned two by two and smoothed with a double pass 3x3 box 
car filter, which eliminates sharp peaks and cosmic rays and minimizes ccntroid errors. Centroids, with sub-pixel 
accuracy, were obtained by a linear interpolation around the brightest pixel in each sub-aperture. The results 
from this fast and simple centroiding method differed from a more precise, moment centroiding analysis by less 
than 0.08 pixels rms. A reference image was built from roughly six thousand frames (60 seconds of data) and 
used to define the zero-deviation, reference positions for the Shack-Hartmann analysis. 

A 36-mode Zernike decomposition (complete up to 6th radial+azimuthal order) was performed on each 
instantaneous image from the Shack-Hartmann imaging. The end product of the data reduction is a time-series 
of Zernike coefficients for each camera. Typically 12,000 frames (2 minutes) of data are analyzed. Note that 
the tip/tilt terms include both atmospheric effects and tracking errors of the mount, particularly wind-shake. 
During the observations guide corrections were sent to the mount using the standard guide cameras in the usual 
way, with a response time of about one second. 

Since GLAO reduces the wavefront error but typically does not result in diffraction-limited images, it is useful 
to be able to consider the effect of the Zernike modes as quadrature contributions to the RMS image diameter 
in arcseconds. Higher-order modes produce larger images for a given wavefront error because the slope errors 
occur over shorter scales. Table presents a list of conversion factors expressed in arcseconds per micron of 
wavefront error. The table was derived using Zernike phase screens in a Zemax model of a perfect 6.5 meter 
telescope. Note that the RMS image diameter is determined from the image centroid, which in the case of coma 
and higher-order coma may neglect some of the error which shows up as a shift in the image centroid. 

Table 3 is a summary of the properties of the Zernike analysis of each run. Column 2 is the seeing as observed 
through the 59 cm sub- apertures. This number is derived from the total wavefront error (col 3), the fractional 
wavefront errors observed in each Zernike order (col 4-10) and the conversion factors given in Table 2. 

Figure n shows an example of the RMS of the Zernike amplitudes for Run 11. The overall shape of the RMS 
as a function of mode is predicted by atmospheric theory and is useful for deriving measurements of r D , the 
outer scale, and testing the validity of Kolmogorov turbulence (to be done in a future study in combination with 
full presentation of seeing monitor data). Note that there is some variation in the overall RMS from camera 
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Table 2. Image Diameter RMS (arcsec) per Micron of Wavefront Error. 
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Figure 1. Example of the RMS of the Zernike modes in nm of wavefront error. The four different symbols represent the 
four cameras from Run 11. 

to camera. This is due partly to fainter stars having noisier measurements, but also partly to small variations 
in the spacing of the optics in the wavefront sensors. The optical effect also includes focus errors caused by 
the curvature of the telescope focal plane, although we attempt to correct for this during the setup for each 
observation. To move the measurements on a common scale, the values for each wavefront sensor are scaled to 
the average for all wavefront sensors in a data set. 

The RMS as a function of Zernike mode is plotted in arcseconds in Figure [3 Most of the modes contribute 
comparably to the total seeing, although it is clear that the contribution is falling off after mode 30 and the 
radial terms (4, 11, 22) contribute more than the average mode. Also note that although 75% of the wavefront 
variance is contained in tip/tilt (Table 3), the contribution to RMS image size from these two modes is quite 
small. 
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Table 3. Total Wavefront Aberrations and Modal Contribution to Total. 
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Figure 2. Example of the contribution to the total image diameter of the Zernike modes in arcsec. The four different 
symbols represent the four cameras from Run 11. 

Figure|21shows the seeing from the Clay guide camera compared with the seeing determined from the wavefront 
sensors. On average, the 59 cm sub-apertures of our wavefront sensors detect 86% of the variance of the total 
seeing. This implies that any fractional ground-layer contribution that we measure applies only to the 86% of 
the seeing that we detect. Presumably the "missing seeing" is in higher order modes, which will have small 
isoplanatic scale lengths and little impact in a GLAO system. 

5. GLAO ANALYSIS 

For each pair of cameras, we decompose the wavefront errors in each Zernike mode into a correlated component, 
Zg-i which we attribute to the ground-layer and an uncorrelated component, Zjj, which we attribute to the 
upper atmosphere. Any noise associated with the system is also attributed to Zu, since the upper atmosphere 
is effectively noise for a system attempting to sense the ground-layer. Then for wavefront sensors 1 and 2, the 
wavefront errors are Z\ = Zq + Zu,i and Z2 = Zq + Zu,2, where Zu,# is the upper atmosphere turbulence 
seen through each wavefront sensor and Zq is common to both cameras by definition. Define P = Z\ + Z2, 
and M = Z\ — Z 2 . The variance of P is < P 2 >= 4 < Z G > +2 < Zfj > and the variance of M is 
< M 2 >= 2 < Z 2 j >. The fraction of the variance arising from the ground-layer is GL = <pii>^<jv/ii> = 

<Zq> Ground— Layer V ariance 

<Z'q> + < Z'fj > ~ Total Variance ' 

In Table 4 we report the fraction of atmospheric wavefront variance seen in the ground-layer as a function of 
Zernike radial order. Where multiple baselines were observed within a run, each is reported as a separate line. 
Tip and tilt are not distinct from mount shake and may not reflect the ground-layer seeing. In the high wind 
data, Runs 13-15, a clear elongation is seen of the spots, indicating some mount shake. In the next section we 
put these eleven nights of ground-layer into context and convert these fractional variance measures into image 
size. As a consistency check, two observations from separate nights were analyzed with this method. These data 
which are known to be completely uncorrelated produced results of GL < 0.02 for all orders. We note that Run 
8 is very close to zero and not consistent with the other runs in any of the modes. We are in the process of 
reexamining this dataset more closely. 
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Figure 3. The Clay seeing is plotted against the seeing detected in the 59 cm sub-apertures of the GLAO experiment. 
The dotted line represents the 1 1 relation. All seeing measurements have been corrected to a reference wavelength of 
0.5/im. The 59 cm sub-apertures detect 86% of the total atmospheric variance. 




Figure 4. Example of the field effect of the ground-layer. The six baselines are plotted for Run 11. Modes are averaged 
together to provide increased signal-to-noise. 

In Figure 01 the fraction of variance in the ground-layer, GL, is plotted against separation for Run 11. The 
expected behavior is that the ground-layer contribution will fall off as field angle to the 5/3 power. However 
the ground-layer fraction appears to fall off more slowly or even remain approximately constant as field angle 
increases. 
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0.54 


0.44 


0.32 


0.18 


0.519 
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0.69 
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0.54 
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Table 4. Fraction of Variance Detected in Ground-Layer 
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Figure 5. Cumulative seeing histogram for Clay (solid), DIMM (dash), MASS (dot-dash) for one- year's worth of data. 
See Table 1^1 for quartile seeing values. 

6. SUPPORTING DATA 

In order to put the GLAO experiment data into context, we have examined a year's worth of site monitoring 
data. Meteorological data are being collected at the Magellan site using weather stations manufactured by Davis 
Instruments Corp. A Weather Monitor II station has been in continuous operation at site since late-2000. In 
2005 this station was replaced with two newer models (Vantage Pro & Vantage Pro2) which utilize the same 
basic sensor technology to measure temperature, humidity, atmospheric pressure, and wind speed and direction. 
The stations generate measurements for temperature, pressure, humidity, wind speed, and wind direction every 
1-2 seconds and are stored in a database as 1 minute averages. 

The seeing is being monitored at the Magellan telescopes through the use of a differential image motion 
monitor (DIMM). The DIMM 11 functions by using Kolmogorov turbulence theory to relate the FWHM of a 
long exposure in a large telescope to the variance in the difference in the positions of two images of the same 
star created by placing two weak prisms at the aperture of a small telescope. The Magellan DIMM is based on 
the CTIO RoboDIMM, but has several improvements. Following a technique developed by the TMT project, 
software was written to use a SBIG ST7 CCD in a drift scan readout mode which allows for many more image 
motion measurements to be made per minute and thus improved statistics. Image quality has also been improved 
by using two thinner prisms as opposed to one thicker prism and an open aperture. 

Turbulence in the atmosphere above 500 m is being monitored by a multi-aperture scintillation sensor 
(MASS, 12 ). The spatial scale of the scintillation variation depends on the distance to the turbulence which 
gives rise to the wavefront phase disturbance. The turbulence distribution in a small number of discrete layers 
can be determined by fitting a model to the differences between the scintillation indices within four concentric 
apertures. Our MASS has an accompanying DIMM built into it. This instrument, known as a MASS-DIMM, 
was fabricated and provided by CTIO and put into operation in a tower at the Magellan telescopes site. The 
MASS-DIMM instrument operates at the back-end of a 0.25 meter Meade telescope in its own, ventilated IO- 
meter tower. The MASS is only sensitive to turbulence above 500 m. In addition to the turbulence profile, 
the MASS also measures free atmospheric seeing (essentially the integral of the turbulence profile) , the adaptive 
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0.84 


10% 


1.07 


1.07 


1.07 



Table 5. Seeing Percentiles for MASS, DIMM, and Clay 



optics time constant and the isoplanatic angle. The difference between the DIMM seeing and the free atmosphere 
MASS seeing is a measure of the seeing contributed by the ground layer. 

In addition to the site monitoring equipment, the seeing from both Magellan Telescopes is recorded auto- 
matically. The Magellan guide cameras have low quality optics with significant aberrations and so a correction 
(0.35" in quadrature) must be removed from the observed FWHM values. The guide cameras use a RG610 filter 
and that combined with the midband coated E2V CCD47-20 produces an effective wavelength of 765 nm. 

We examined 205 nights from March 2004 to April 2005 that had both seeing measurements from all the 
instruments and meteorological data. (Note that the clear fraction of nights is higher, but the MASS-DIMM 
instrument was being commissioned in the first half of 2004.) For each seeing measurement, a mean tempera- 
ture, pressure, wind speed and wind direction was defined as the mean meteorological condition for 15 minutes 
surrounding the seeing measurement. In the end, we assembled 64,000 MASS and 64,000 DIMM seeing mea- 
surements and 128,000 Clay guide camera seeing measurements. All of the data were corrected to an airmass of 
1, and to a reference wavelength of 0.5/im. The MASS data were reprocessed with the most recent version of 
the turbina software. 

The first consistency check is to compare the Magellan seeing measurements to the DIMM seeing. In Figure 
El the cumulative seeing histograms are plotted for the Clay guide camera (solid), DIMM (dash), and MASS 
(dot-dash) and in Table HO select values are presented numerically. The agreement between the DIMM and Clay 
is good, but not perfect. The low end difference between the two is likely caused by the simple correction applied 
to compensate for the guide camera optics. This correction is a first order correction and could likely be improved 
with further analysis. The effects of an outer scale will suppress the seeing in the larger Clay Telescope aperture. 
However, since the DIMM is a differential measure (and a much smaller structure), the DIMM seeing values will 
be unaffected by mount and wind shake. It appears that to first order, these two effects cancel each other out 
and we obtain good agreement between the Magellan guiders and the DIMM. 

The MASS data relative to the DIMM data indicates that in the median conditions, ground-layer is responsible 
for about 38% of the atmospheric variance. This is a smaller percentage of total seeing than is observed at other 
sites. 4 ' 13-15 At Las Campanas, a ground-layer correction would work in the sense of making a good night better 
as there is little difference seen between the MASS and DIMM at higher seeing values. 

In Figure El the ground-layer as measured by MASS-DIMM, (1 — jy£r§p ) , is plotted against the trends lines 
for ground-layer fractional variance as measured by our experiment. We set any MASS-DIMM data that results 
in a negative value to zero and attribute this to noise and small errors in the absolute calibration between the 
two instruments. There is a correlation in between the MASS-DIMM and the GLAO experiment, although the 
relationship is far from one-to-one. The crosses are the data that were used to determine the trend-line for orders 
3-6 and provide a sense of the uncertainties. It is interesting to note that even when the MASS-DIMM indicates 
that there is no ground-layer, there is a significant tip/tilt correlation, indicating that approximately 30% of the 
variance in tip/tilt is due to mount shake. 

There is also a small offset in the higher order modes. Since a GL measurement between two separate nights 
produced a result of less than 0.02, the positive values of about 0.1 in the higher-order modes, even when the 
MASS-DIMM comparison fails to detect the ground layer, are likely to be significant. 

Figure [3 shows the ground-layer as measured by MASS-DIMM and the RMS image improvement. These 
numbers were derived by reducing each mode in the RSS of the total variance by the amount of observed 



Ground Layer Fractional Variance from MASS-DIMM 



Figure 6. MASS-DIMM measured ground-layer vs ground-layer from this experiment. The mean relation between 
ground-layer variance as measured by the MASS-DIMM instrument and the Magellan Telescope is plotted as a solid line 
and labeled All Modes. The data for All Modes are plotted as crosses to provide a sense of the errors. Trend-lines for 
tip/tilt (dot), Focus, Astigmatism, Coma, Trifoil (dashed), and higher orders (dot-dash) are drawn as a fit to the data 
from Table |5] The vertical lines indicate the quartile conditions as determined from a year's worth of seeing data. 

ground-layer reduction (Table 4). These scaled variances are then converted to arcseconds and compared to the 
uncorrected RSS variance from the same data. These calculations include the scaling factor that accounts for 
the fact that the 59 cm sub-apertures only see 86% of the seeing. 

The largest gains are obtained when the MASS-DIMM reports the strongest ground-layer. Note that multiple 
points will show up for a given observation if more than one baseline was observed. It is clear that strong ground- 
layer nights allow for improvements in image quality, although the amount of improvement is significantly less 
than predicted. Neither of the two previous figures are changed significantly by the removal of the largest (>10 
arcmin) baseline data. The strongest ground-layer night also happened to be the night of the strongest winds 
(and technically over the operating limit!). However, no correlation exists between image correction and ground 
wind speed. 

7. CONCLUSION 

We have taken one of the first sets of moderate speed (100 Hz), wide-field (arcminutes), and moderate resolution 
(59 cm sub-aperture) wavefront sensor data, in order to asses the gains of a GLAO system. We have taken 
this data on eleven nights which span a range of meteorological, seeing, and reported ground-layer strength 
conditions. On the best nights, we find a RMS image size reduction of 30% over a 7 arcminute separation. On 
a typical night we find a 10% correction to a visible-band image. 
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Figure 7. Improvement in RMS images diameter shown as a function of MASS-DIMM reported ground-layer. Vertical 
lines indicate the quartile and median conditions. Note that not all the runs have MASS-DIMM data, and each baseline 
for a given run is plotted, providing multiple measurements for some nights. This calculation includes a correction to 
account for the fact that 86% of the total seeing variance is observed in the 59 cm sub-apertures. 
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